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Abstract 

Oscillations of three pseudo-Dirac flavor neutrinos z/g, z/^, z/^ are considered: < 
^(i) _ ^(R) ^ m^^^ for their Majorana and Dirac masses taken as universal before 
family mixing. The actual neutrino mass matrix is assumed to be the tensor product 
M^"^ (g) I ^ ^(j^) j, where M^"^ is a neutrino family mass matrix (M^'')"^ = M^"")) 
and \^^'^^ = m^^'^ym^^\ The M^'^^ is tried in a form proposed previously for charged 
leptons e, fi, r for which it gives m,- = 1776.80 MeV versus m°^P = 1777.051^020 MeV 
(with the experimental values of and used as inputs). However, in contrast to the 
charged-lepton case, in the neutrino case its off-diagonal entries dominate over diagonal. 
Then, it is shown that three neutrino effects (the deflcits of solar z/g's and atmospheric 
z/^'s as well as the possible LSND excess of u^s in accelerator z/^ beam) can be explained 
by neutrino oscillations though, alternatively, the LSND effect may be eliminated (by a 
parameter choice). Atmospheric z/^'s oscillate dominantly into z/^-'s, while solar z/g's — 
into (automatically existing) Majorana sterile counterparts of u^s. 
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Let us consider three flavor neutrinos z/g, ^'^t, and assume for them the mass matrix 
in the form of tensor product of the neutrino family 3x3 mass matrix (jXlj^^^ [a,f3 ~ 
e, II, r) and the Majorana 2x2 mass matrix 

the second divided by m^^^ (with m^^^ included into m'^^). Then, the neutrino mass 
term in the lagrangian gets the form 



1 /-W -W\ (,) 




where 



o(a) o /o \c o(s) o /o \c 

l^a =^aL + [l'aL) , I'a = R + [^a r) (3) 

and A^^'''^^ = m^^'^^ /m^^\ Here, and I'U are the conventional Majorana active and 
sterile neutrinos of three families as they appear in the lagrangian before diagonalization of 
neutrino and charged-lepton family mass matrices. Due to the relation v^^/s — ly^i^a, the 
family mass matrix M^''^ = M^^^^^, when standing at the position of A^-'^) and A^^^ in Eq. 
(2), reduces to its symmetric part |(M('')+M('')^) equal to its real part KM^^^+M^*") *) = 
ReM'^'^K We will simply assume that (at least approximately) M'-"^ = M^'')'^ = M^'')*, 
and hence [/^ = [/(^)* = (t/M-^) . Then, CP violation for neutrinos does not appear 
if, in addition, C/(^) = C/(^)*. Further on, we will always assume that < A^^^ = A^^^ (= 
A(^)) and A^^) < 1 (the pseudo-Dirac case) [1]. 

Then, diagonalizing the neutrino mass matrix, we obtain from Eq. (2) 

->c_ = ^EH, ( A- ) ( ' 

where 
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(i, j = 1,2,3) and 

with Via — {u^'^^ , ^ Up'l describing the lepton counterpart of the Cabibbo — Kobayashi — 
Maskawa matrix. Here, 

and 

where (m^^^) (q;,/9 = e, r) is the mass matrix for three charged leptons e~, //~, r~, 
giving their masses me, m^, after its diagonalization is carried out. Now, u^f' and u'^^ 
are the conventional Majorana active and sterile flavor neutrinos of three families, while 
vl and are Majorana massive neutrinos. 

If CP violation for neutrinos does not appear or can be neglected, the probabilities 
for oscillations — >■ i/^"'* and u'^^ — >■ u'^i^^ are given by the following formulae (in the 
pseudo-Dirac case): 



j>i 



sin^ (xj-xl) + sin^ (xf -x-^) + sin^ (xf -x-) + sin^ 



(9) 



and 



P ( i^i"^ 



-EVjpVUVrpVia [sin^ (xj-x[)+sin2 (xf -xf ) -sin^ (xf -sin^ (^-^f ) 



(10) 
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where P|i//'-^^) = pf'^V/'-^^) , p^' = ^E^ - (m,,A-f.^^)2 E - {m,,X''''f/2E and 



1^2 \I,II\2t 

1,11 ^ 1.27 ^ '-'^ ^ , = 1 T 2AW ^ 1 (11) 

E 



with mj/., L and E expressed in eV, km and GeV, respectively (L is the experimental 
baseline). Here, due to Eqs. (11), 



-x\ = \.n^^ — (12) 



and for j > i 



Then, the bracket [ ] in Eq. (9) and (10) is reduced to 4sin^ \.TI{m\. - ml.)L/E and 0, 
respectively. The probability sum rule X)/3 P {ya^ ~^ ^^P^^ ^ (^a^ ~^ ^/9*0 ~ ''" 
readily from Eqs. (9) and (10). 

Notice that in the case of lepton Cabibbo — Kobayashi — Maskawa matrix being nearly 
unit, (Via) — (Sia), the oscillations — > v^^^ and — > u^^^ are essentially described 
by the formulae 



(a) 
/3 



J /3a 



P(v^:^^vf) ^ 5,„sin^(l.27^ 



4ml A^L^ 



E 



(14) 



corresponding to three maximal mixings of with i/^^ {a — e, r). Of course, for a 
further discussion of the oscillation formulae (9) and (10), in particular those for appear- 
ance modes z/^"-* v^^ (a 7^ /?), a detailed knowledge of (V^q,) is necessary. 

To this end, we will try to extend to neutrinos the form of charged-lepton mass matrix 



2Q;(^)e-^^^'' 4/x(")(80 + £(^))/9 8^3 Q;(^)e^^''^ 
\ Sv^a^^^e-^*'*'^ 24//(^)(624 + £(^))/25 / 



(15) 
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which reproduces surprisingly well the charged-lepton masses mg, m^, o;^^^ and 

s^^) are positive parameters). In fact, treating off-diagonal elements of (m^^''^ as a per- 
turbation of its diagonal entries, we get the mass sum rule 



6 / a'^"^ 

(351m„ - 136me) + 10.2112 ^ 

125 ^ ' V//(^) 



MeV 



1776.80+ 10.2112 



a 



MeV , 



(16) 



where the experimental values of mg and are used as inputs. Then, /x^^^ — 85.9924 
MeV and s^^^ = 0.172329 (up to the perturbation). The prediction (16) agrees very 
well with the experimental figure = 1777.05lo;|o MeV, even in the zero order in 

. Taking this experimental value of rrir as another input, we obtain 



a 



n24+° °2^ 



(17) 



what is not inconsistent with zero. 

Now, we conjecture the neutrino family mass matrix (^M^^^^ in the form (15) with 

s^'^) ~ and ip^^^ cp^''^ = [2]. In order to get the 



neutrino family diagonalizing matrix [uj;^-^ rather different from the unit matrix (S^^} j, 
we assume that diagonal elements of (^M^^^ can be considered as a perturbation of its off- 
diagonal entries (though the diagonal as well as the off-diagonal elements are expected 
to be very small). Under this assumption we derive the unitary matrix (Uai^ of the 
following form : 



with (p^"^ = and 



/ Vis 

7 



1 



-2i(^M Vis 



1 

V2 



4S „-i(p'-''^ 



7V2^ 



73' 



A/is 
7V2 



+ 0{C/7) 



(18) 



M 



33 



299.52 



a 



X 



16.848 



(19) 
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In this case, the neutrino family masses are 



m 



1^1 



49 1 12 I 



T7 + 



1 /48 



2 V49 



\Mii>\ , 



(20) 



where \Mii>\ = 2q;('^)/29 (thus, m,^i <^ |m,^2l < "^i^s)- Hence, 



< - < = 14 + x) |m£) p = 20.721 



(21) 



Taking in contrast (Uap) — i^a/i) — as in (^M^^^ the off-diagonal elements are pertur- 
batively small versus diagonal entries [cf. Eq. (17)] — we can insert 



(u)* 



(22) 



into Eqs. (9) and (10). Here, U^i are determined from Eq. (18) 



Then, with the use of Eqs. (12) and (13) the oscillation formulae (9) take 

the form 



p ^^ia) ^ ^(a)) 



l_llsin2 (1.27- 



492 
1 



E 



4-492 
96 



2 , , _4m^,AWL\ 



sin2 1.27-^ 



E 



492 
1 



sin^ 1.27^^ ''^^ 



■sin^ 1.27- 



492"" ■ ^ 
sin^ 1.27 ^ 



+ sin^ 1.27- 



4miAWL\ 



^-3 



+ sin^ i i.27^^^2--^^^ 



1-1 

4 



4m2 A(^)L\ 
+ sin^ ( 1.27 J 



sin" I 1.27^^ 
1 



E 



4-49 



sin2 1.27—^ 



E 



+ sin^ 1.27- 



4miAWL^ 



1^3 



E 



+— sin^ 1.27^^^2 ^ 

49 I E 



(23) 



and 
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48 



4-49 

48 



sin^ 1.27 j + sin^ ( 1.27 

2 f. n^i^l3-^l2)L\ 



+ — sinMl.27 
49 \^ E 

sm^ 1.27 ^ I 



492 

4; 



4-492 

96 
492 
48 



sin^ 1.27- 



E J 

Ami A(-'^')I~ 



E 



+ 



,2(1 O'T ^ ^^2 ^'1'' 



sin^ 1.27- 



492 



sin^ 1.27- 



E 



-h sin^ 1.27- 



+ sin^ 1.27- 



4//;2 A(^^)L^ 



E 

ml,-ml)L 



E 



E 



1 „/ 4m2AWL^ 
1 - — sm2 1.27- ^1 



492 

482 



4-492 
96 



sin2 1.27^^ 



E J 



492 
482 
"492 



sin^ 1.27- 



E 



2 / 4m2 aWl\ 
+ sin^ [1.27^ j 

iml^-mljL' 



E 



+ sin2 1.27- 



E 



(24) 



In these formulae, the experimental baselines L (and neutrino energies E) are generally 
different. 

Further on, wc intend to relate the first, second and third Eq. (23) to the experimental 
results concerning the deficit of solar i/g's [3], the deficit of atmospheric i/^'s [4] and the 
excess of i/gS in accelerator beam [5], respectively. 

First, let us assume the simphfying hypothesis that the LSND effect [5] does not exist. 
Then, under the numerical conjecture that 



^o /<AWL,,i ^^^^^ ^ ^3 /<AWL,,^ ^ ^/ m2^W-gatm \ 

Esoi Ea,tm \ rnl^Lsoi/ Esol J 

1.27^<~<^^"-^0(1) , 1.27^^-3 -O^soi ^ of^^^4^')»l, (25) 
we obtain from Eqs. (23) 
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49^"" I'-'' ^soi J - 4:49^ 
482 . 2^^ ^^4m^^AWL,oi\ 



1-— sin^ 1.27 

492 I K 



sol 



atm / -'^atm 

2 ^2 



2 / 1 „„(?7T.j,3 — 771^,2 )-Z^atm 



~ l-sin^ 1.27^^^^^^^ ^ 



^atm 



P(,0.)^,(a)) ^ -^(1.27)2 -"--% "^^^^ +^(1.27)^ 



49^--^ ^EsND ■ 49^— Els^^ 
~ . (26) 

The term —387/4 • 49^ = —0.0403 in the first Eq. (26) comes out from averaging all sin^ 
of large phases over oscillation lengths defined by sin^ of a phase = 0(1) (then, each sin^ 
of a large phase gives 1/2). 

Comparing Eqs. (26) with experimental estimates, we get for solar i^e's [3] (using the 
global vacuum fit) 

48^ 

— ^ sin^ 2^,01 - 0.75 , Aml^X^^^ ^ Am^^i - 6.5 x 10"^^ eV^ , (27) 
49 

and for atmospheric i/^'s [4] 



1 ^ sin^ 2^atm - 1 , ml^- ml^ ^ Am^,^ - 2.2 x lO'^ eV^ . (28) 



Thus, from Eqs. (27) and (28) 



^.3.0x10-. (29) 



Hence, making use of Eqs. (20) and (21), we infer that 



- 2.6 X 10-^ , Il^xi^) = l^A(^) ^ 7.5 x 10"^ , 

= TTT^eA^"") - 8.6 X 10-^ , a(^V^"^ = ^^7^ 1-1 x 10"^ eV^ , 

299.52^ ' ^ 20.721 

//M2^(^)~ 9.1 X 10-^1 eV2 , ^^i.2xl02eV2 . (30) 
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Here, the constant ^ still may be treated as a free parameter (determining A*^^)). If 
^ = O(lO-i), then = 0(10-^), m^J\m^^\ = 0(10-^), ii^^^la^""^ = 0(10-^), At^*^) = 
0(10"^ eV), a(^) = 0(1 eV) and 

= 0(10-^ eV) , |m^,| = 0(10"^ eV) , m^^ = 0(10"^ eV) (31) 

with ml^ - ml^ ~ 2.2 x lO^^ eV^. 

In this way, both neutrino deficits can be explained by pseudo-Dirac neutrino oscil- 
lations. Note that solar i/^")'s and atmospheric ^'^"^'s oscillate dominantly into i/^^^'s and 
4"^'s, respectively (here, i/^"], = i^aL, i^aL = (^^a)L )■ 

Now, let us accept the LSND effect [5]. Then, making the numerical conjecture that 



-E-sol -E'atm Lso\/ Esol J 

-E^LSND -E'atm \ -f'LSND/-E'LSND / 

(32) 

we get from Eqs. (23) 



P(^(a)^^(a)) ^ l-l|sinMl.27 



''%in- o, 4<A(^)L.oA 387 

49^'^" l'-'^ j 4.492 



. l-l^sin^fl.27-^<^^"^^-^ 



492 V ■ K 



sol 

'M) T 

atm 



2 2 V Eg^iYn 



^ ^ 4y \ -C/T.SND / 4y £/T 



LSND / -^LSND 



49 V £;lsnd y ^ ^ 

When comparing Eqs. (33) with experimental estimates, we obtain for solar UgS [3] 
(making use of global vacuum fit) 

^ ^ sin^ 2^,01 ~ 0.75 , 4m^^ AW ^ Am,'„i ~ 6.5 x 10"^^ eV' , (34) 
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for atmospheric f^'s [4] 



1 ^ sin^ 261; 



atm ~ 1 , 2 + 2 



E.. 



sm 



1 



atm 



(35) 



atm 



with 



Am^,^ ~ 2.2 X 10-=^ eV^ , (36) 

and for accelerator i/^'s [5], say, 



^ ^ sin^ 2^LSND ~ 0.02 , - ml^ ^ Am^sND ~ 0.5 eV' . (37) 



So, from Eqs. (34) and (37) 



4^2 \{.M) A 2 



Hence, due to Eqs. (20) and (21), 



~ 1.1 X 10-^ , p^X^m ^ 3.3 X 10-9 , 
~ 3.8 X lO-'' , a^'^)//^ ~ 2.4 x 10"^ eV^ , 



a 



2a W ^9.2 X 10-^1 eV^ , ^ 6.3 x 10^ eV^ . (39) 

A*^ 

Here, the constant ^ still may play the role of a free parameter (determining A^^^). If 
i = O(lO-i), then AW = 0(10-^), m^J\m^^\ = 0(10"^), //M/aM = 0(10"^), /x^ = 
0(10-3 eV), = 0(10 eV), and hence 

m^i = 0(10-3 eV) , Im^^l = 0(1 eV) , m^g = 0(1 eV) (40) 
with — 0.5 eV^. Then, in Eq. (35) we can put approximately 



2 \ Egxm ) 



06 
05 



(41) 
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in a reasonable consistency with the Super-Kamiokande estimate [4] . Here, {U — D)/{U + 
D) is the up-down assymetry for f^'s, estimated as —0.296 ± 0.048 ± 0.01. 

In this way, therefore, all three neutrino effects can be explained by pscudo-Dirac 
neutrino oscillations. Note that solar i/("^'s and atmospheric i/^"'''s oscillate dominantly 
into f^^^'s and v'f^^s, respectively, as in the previous case when the LSND effect was 
absent. 

The recently improved upper bound on the effective mass {my^) of the Majorana u'f^ 
neutrino extracted from neutrinoless double P decay experiments is 0.2 eV [6]. In our 
pseudo-Dirac case, this mass is given by the formula (if Via — U^i*)'- 



as 



49 -29 1 49^ ^^y* 473.96 ' 

c^^'^) = in Eq. (18) (here, U^i = U*^ and A^"^^ = ^1 + A^^^, while 



(42) 



-^=i:uS^{^l + -l') ■ (43) 
Thus, in the option excluding or accepting LSND effect we estimate from Eq. (30) or (39) 
that 

, , f 5.4 X lO-^a(^) - 0(10-6 eV) ,,,, 
^"^-^ ~ I 2.4 X lO-^aH ~ O(10-« eV) ' ^^^^ 

respectively. Thus, in this pseudo-Dirac case, the Oi//3/3 decay violating the lepton number 

conservation is negligible. Note that {rrii,^) <^ rrii,^ <^ \mv2\ < "^i^s iii both options. Here, 

the neutrino masses are 



mif = m,\'^'' = m,, (^1 + A^^)) ~ ^m,, . (45) 

Since for relativistic particles only masses squared are relevant, the " phenomenological" 
neutrino masses are equal to ~ |mi,. | i.e., ~ m^-^ , Im^^l , m^^. 
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Finally, let us turn back to the option, where there is no LSND effect. In this case, the 
natural possibility seems to be a (nearly) diagonal form of neutrino family mass matrix 

jVf('^) ~ {5ai3mu^) and so, unit neutrino diagonalizing matrix U'^^'' ~ {5ai)- Then, if 
U^^^ ~ {5^p) i.e., V ~ {5i^), Eqs. (14) hold, giving 



P {ut^ - ui^^) 1-P ^ .W) c 1 - sin^ (l.27^<^^ , 

1.27^ 



(46) 



Here, rriui = rriy^ are neutrino family masses. 

Comparing Eqs. (46) with experimental estimates for solar i/g's [3] (using the global 
vacuum fit) and atmospheric f^'s [4], we have Eq. (27) (with m^,. = fn^^) and the relation 

1 ^ sin^ 2^atm ~ 1 , ^ml\^^^ ^ Aml^ ~ 2.2 x IQ-^ eV' , (47) 
respectively. Hence, 

^ ^ 3.0 X lO-s . (48) 

Under the conjecture that M^"^ has the form (15) with pi^^^ //(^), a^^) a^"^ = 0, 
— > g^*^) ~ 0, we get 

^(i^ ~1.7xl0-^ (49) 



mj. 4 • 80 yAm^ 



atm , 



Then, 



- 6.1 X 10"^ (50) 



and 



m.^ = ^/x(^) ~ 2.1 X lO-V^"^^ , 
m., = ^ /x^'^) = 1.2261 /i^'^) , m,, = /.^ = 20.657 //^ = 16.848m,^ . 



(51) 
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Here, the neutrino masses are ml'J^ — m,^^A^'^^ = m,^^(^l + A^^^) ~ l^nii,^, so that 
\ml'J^\ ~ m^^ , rriy^ , rriy^ , where m^^ : m^^ : m,^^ ~ 1.7 x 10'^ : 1 : 16.8. Prom Eqs. 
(47) and (51) we infer that 

;,(-)2aW ^3.7x lO-^eV^ . (52) 

In this way, both neutrino deficits can be explained by osciUations of unmixed pscudo- 
Dirac neutrinos {U^^^ ~ 5ai)- Note, however, that now both solar i/(")'s and atmospheric 
i/^"^'s oscillate dominantly into Majorana sterile neutrinos: f^^^'s and J^^^^'s, respectively 
(in contrast to the previous mixed pseudo-Dirac u^"'^ and uj^^ neutrinos of which the latter 
oscillated dominantly into uif'^^s). The experimental evidence for u^- oscillations 

and/or for the LSND effect would be, of course, crucial in the process of understanding 
the mechanism of neutrino oscillations. 

In the present case, the effective mass {m^^} of the Majorana i^^"^ neutrino is given as 

(m.,) ^ m.,^ (y + A^^) = m^^AW , (53) 

since ~ 6ai- Thus, the Ovfifi decay upper bound {niy^) < 0.2 eV is certainly satisfied 
because of A^'^^ <^ 1 (and rrii,^ < a few eV). 

If it turned out that both solar i/g's and atmospheric f^'s oscillated into sterile neu- 
trinos, it would not be easy to recognize whether, as discussed above, the latter should 
be Majorana sterile counterparts of Majorana active i^gS and I'n^s, or rather, two extra 
Dirac sterile neutrinos [7]. 
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